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Abstract 

Dark matter may be hidden, with no standard model gauge interactions. At the same time, in 
WIMPless models with hidden matter masses proportional to hidden gauge couplings squared, the 
hidden dark matter's thermal relic density may naturally be in the right range, preserving the key 
quantitative virtue of WIMPs. We consider this possibility in detail. We first determine model- 
independent constraints on hidden sectors from Big Bang nucleosynthesis and the cosmic microwave 
background. Contrary to conventional wisdom, large hidden sectors are easily accommodated. A 
flavour-free version of the standard model is allowed if the hidden sector is just 30% colder than 
the observable sector after reheating. Alternatively, if the hidden sector contains a 1-generation 
version of the standard model with characteristic mass scale below 1 MeV, even identical reheating 
temperatures are allowed. We then analyze hidden sector freezeout in detail for a concrete model, 
solving the Boltzmann equation numerically and understanding the results from both observable 
and hidden sector points of view. We find that WIMPless dark matter indeed obtains the correct 
relic density for masses in the range keV ^ mx TeV. The upper bound results from the 
requirement of perturbativity, and the lower bound assumes that the observable and hidden sectors 
reheat to the same temperature and is raised to the MeV scale if the hidden sector is 10 times 
colder. WIMPless dark matter therefore generalizes the WIMP paradigm to the largest mass range 
possible for viable thermal relics and provides a unified framework for exploring dark matter signals 
across nine orders of magnitude in dark matter mass. 

PACS numbers: 95.35.+d, 04.65. +e, 12.60.Jv 
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I. INTRODUCTION 



At present, all solid evidence for dark matter is gravitational. At the same time, the possi- 
bility that dark matter has electromagnetic or strong interactions is highly constrained [l|, 0] ■ 
A straightforward possibility, then, is that dark matter is hidden, consisting of particles that 
have no standard model (SM) gauge interactions. 

Hidden sectors have a long and distinguished history. For example, the idea that a hidden 
sector may restore parity on cosmological scales is as old as the idea of parity violation 
itself [3||. Such hidden sectors, containing "mirror matter," have been studied by many 
groups, as have their possible implications for dark matter [1, M, @, 0, H, @, EJ 11 Jl2l 13 . 
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String theory also motivates hidden sectors with their own gauge interactions 14| 
the case of the heterotic string [l5[ and intersecting brane models [l6|. More recently, 
hidden sectors have been central to several phenomenological developments, including Higgs 
portals [l7, 18], hidden valleys [l9[, unparticles [13], and quirks [2lj], and the possibility that 
dark matter may reside in a hidden sector has been discussed in several recent works [22I, 

2!, 0, B M, m H H H 51. 



In considering hidden dark matter, one would seemingly be sacrificing critical virtues of 
more conventional dark matter candidates, namely, the connection of dark matter to the 
gauge hierarchy problem, and the fact that weakly-interacting massive particles (WIMPs) 
naturally have the desired thermal relic density. This is not necessarily true, however. In the 
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recently proposed framework of WIMPless dark matter [30| , hidden sector dark matter also 
naturally has the desired thermal relic density. The essential idea is that, very generally, 
a dark matter candidate that decoupled from the thermal bath non-relativistically has a 
thermal relic density 

1 m 2 

where (ctav) is the thermally-averaged annihilation cross section times velocity, and g and m 
are the dark matter particle's gauge coupling and mass. For WIMPs with weak interaction 
coupling constant g wea k — 0.65 and weak-scale mass m wea k ~ 100 GeV — 1 TeV, this relic 
mass density is naturally near the observed Odm^ 2 — 0.11. 

In WIMPless models, hidden sector particles X have masses and couplings that may be 
very different from WIMPs, but which nevertheless satisfy 
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In the WIMPless examples discussed in Ref. [30|, Eq. ([2]) follows from the structure of 
gauge-mediated supersymmetry breaking (GMSB), which generates hidden sector masses 
that are proportional to couplings squared. The WIMPless framework therefore naturally 
generalizes the WIMP paradigm to other cold dark matter candidates without sacrificing the 
key thermal relic density virtue of WIMPs. If additional connector particles with both SM 
and hidden gauge quantum numbers are present, they will mediate SM-hidden interactions. 
The WIMPless framework may therefore also predict qualitatively new signals, such as an 
explanation of the DAMA signal in terms of GeV dark matter |32|. 

In this paper, we study hidden thermal relics in general, and WIMPless dark matter in 
particular. For the most part, we assume that there are no connector particles, so that the 
hidden sector interacts with the SM only very weakly. In Sec. [Til we carry out a model- 
independent analysis of constraints on hidden sectors from Big Bang nucleosynthesis (BBN) , 
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the cosmic microwave background (CMB), and other cosmological data. As is well known, 



these exclude the possibility that the hidden sector is an exact copy of the SM [33j. This 
statement assumes that the hidden and observable sectors are at the same temperature at 
late times, however. In the present context, there is no strong reason to assume that the 
observable and hidden sectors are reheated to the same temperature (1, @], and, even if they 
are, for a truly decoupled hidden sector, their temperatures will differ by the time of BBN. 
Following previous works in the context of mirror matter (see, e.g., Refs. @, Hi), 
we consider the cosmological constraints allowing different temperatures and find that even 
large hidden sectors are easily accommodated. For example, a hidden version of the SM is 
allowed even if the hidden sector's reheating temperature is as much as 70% of the observable 
sector's. Alternatively, if the hidden sector contains a 1-generation version of the MSSM with 
characteristic mass scale below 1 MeV, even identical reheating temperatures are allowed. 

We then present a concrete model for the hidden sector in Sec. [Till and study the freezeout 
of dark matter in this hidden sector in Sec. HVl The relations of Eqs. ([T]) and §2§ are, of 
course, only rough relations capturing the parametric dependences. In these Sections, we 
present a detailed analysis of a specific case that highlights several subtleties and generalizes 
previous freezeout analyses to cases with sectors at different temperatures. This analysis 
confirms the validity of conclusions based on Eqs. ([I]) and ([2]) and highlights various subtleties 
of the WIMPless framework. In particular, we find that WIMPless dark matter may have 
the desired thermal relic density for the entire range of keV ^ nix ~ TeV, where the lower 
bound is set by the requirement that dark matter freezeout is non-relativistic, and the upper 
bound follows from the requirement of perturbative gauge couplings. WIMPless dark matter 
therefore encompasses as large a range of masses as one could expect of dark matter that has 
the naturally correct thermal relic density, and it provides a unified framework for addressing 
many diverse dark matter signals and phenomenology. 

Finally, in Sec. |V], we discuss modifications from the presence of connector particles 
with both hidden and observable sector gauge interactions. We find that, under general 
assumptions, these fields do not upset the conclusions derived earlier. We present our 
conclusions in Sec. ED 



II. BBN AND CMB CONSTRAINTS 

We assume that the observable sector is the minimal supersymmetric standard model 
(MSSM), which is supplemented by a single hidden sector. Hidden sector parameters are 
denoted by the superscript h. For most of this work, we assume that there are no connectors, 
fields with both SM and hidden gauge charges. The observable and hidden sectors therefore 
interact extremely weakly with each other, and the observable sector's temperature T need 
not equal the hidden sector's temperature T h . We define 

«*) = 5f (3) 

to parametrize the mismatch. In addition, it will be convenient to introduce the standard 
relativistic degrees of freedom parameters. Letting 

_ f 1) i — boson . v 
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we define 



9*(T) = E C i9i (^V 9*s(T) = E C i9i (^V , (5) 
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where the sum is over observable sector particles, and 
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where the sum is over hidden sector particles. Tj is the temperature of particle i, and ^ 
denotes its internal degrees of freedom. Throughout this work, we assume that the gravitino 
has a negligible relic density, as is the case, for example, for low and moderate reheating 
temperatures, and the gravitino is not included Eqs. ([5]) and §6§. 

BBN is sensitive to the expansion rate of the Universe at time £bbn ~ 1 s and temperature 
^bbn ~ 1 MeV. The expansion rate is determined by the energy density, and so constrains 



light particles even if they have no SM interactions [35j, [36|. The constraint is conventionally 
quoted as a bound on iV eff , the effective number of light neutrino species, and may be taken 
to be iV e fj = 3.24 ±1.2 (95% CL), where the baryon density has been fixed to the value 



determined by the CMB, and both He and D data are included [37|, [38j. In our framework, 
this implies 

^(7hb N ) ( ^1 = l ■ 2 ■ (iV eff - 3) < 2.52 (95% CL) , (7) 



BBN 



where Tg BN is the temperature of the hidden sector at time £bbn- If the hidden sector is an 
exact copy of the MSSM and Tg BN = T B bn, 9* (^bbn) = 10.75. This violates the bound of 
Eq. ([7]), as is well-known j33[, and may create the impression that large hidden sectors of 
the size of the MSSM are completely excluded. 

If the observable and hidden sectors are not in thermal contact, however, the hidden 
sector may be colder than the observable sector. This would be the case if, for example, the 
inflaton couplings to the observable and hidden sectors are not identical, so that they reheat 
to different temperatures 0,0]. Alternatively, the observable and hidden sectors may initially 
have the same temperature, either because they have the same inflaton couplings or because 
they are in thermal contact, but may cool independently and have different temperatures 
at later times. In fact, it quite reasonable that the observable and hidden sectors lose 
thermal contact before BBN. Thermal contact would require efficient interactions between 
the two sectors at low temperatures. These would require light connector fields to mediate 
these interactions, but light fields with SM gauge charges generically violate experimental 
constraints. We conclude, then, that it is far from guaranteed that the observable and 
hidden sectors will be at the same temperature at £bbn; generically, they will be at different 
temperatures. 

To explore this possibility, we assume that the observable and hidden sectors reheat at 
time £rh to temperatures T R h and T^ n , respectively. We further assume that the entropy 
per comoving volume is conserved in each sector separately, so that 



g*ff(^BW)^BBN _ g*s(?BBN)?BBN 
9*s ( t rh ) Trh 9*S (Trr ) T| 



(8) 



RH 



Given mass spectra in the observable and hidden sectors, reheat temperatures Trh and T^ H , 
and Eq. (jSJ), the ratio ^bbn /^~bbn is fixed, and one can determine if Eq. ([7]) is satisfied. 
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To provide numerical examples, we define 



rph 

£rh = ^ , (9) 

J- RH 

the ratio of temperatures just after reheating, and assume Trh is above the mass of all 
MSSM particles. We also define 

Slight = X/ CiSi ^heavy = /] CjQj (10) 

for the observable sector and 

^eavT- ^ E Ci* (11) 

for the hidden sector, where we have assumed all hidden relativistic degrees of freedom have 
a common temperature. With these definitions, Eq. (jS]) becomes 



n hBBN /rph \* 

flight / J BBN \ Slight c s 



„/iBBN I „/iBBN I rp I , 

Slight Shcavy \ J BBN / Slight T Shcavy 



e RH • (12) 



Solving for Tg BN /^bbn an d substituting the MSSM values flight = 10.75 and flight + Sheavy 
228.75, we find that the BBN bound of Eq. (ffl) becomes 



/iBBN ^ 3.49 ^ / ft,BBN\l hBBN /io\ 
Sheavy — ^ ^ TT j Vylight J Slight ■ \ LO ) 
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Bounds in the (g^S^^ , Sheaw') pl ane l° r various values of ^rh are given in Fig. [H Several 
aspects deserve comment. As evident from Eq. (1131) . the bound is extremely sensitive to 
£rh- The £rh = 0.5 contour is barely contained on the plot, implying that if the hidden 
sector reheats to half the temperature of the observable sector, very large hidden sectors are 
allowed. Several example hidden sectors, defined precisely in Sec. IHIt are also given on the 
plot. Model A is a 1-generation flavor-free version of the MSSM, with all Yukawa couplings 
of order 1. We see that it is allowed for £ RH < 0.92, that is, reheat temperatures that are 
almost identical to the observable sector. This is our prototype model, which will be used in 
our detailed freezeout analysis in Sec. IIVI Model B is a full 3-generation flavor-free version 
of the MSSM. We see that even a hidden sector as "large" as the MSSM is allowed, provided 
£rh < 0.74. 

Perhaps even more interesting, Fig. [1] shows that Model C, a 1-generation version of 
the MSSM with all degrees of freedom relativistic at BBN, is allowed, even for £rh — 1-06. 
Naively, one might expect a hidden sector that starts at the same temperature as the MSSM 
and has over 100 light degrees of freedom to be completely excluded. In fact, however, such 
a hidden sector is much colder than the MSSM at BBN times, since its cooling is not slowed 
by the disappearance of heavy degrees of freedom, in contrast to the MSSM. This cooling is 
critical, given the 4th power of the temperature in Eq. ([70, and makes such a hidden sector 
allowed. For similar reasons, Fig. [T] shows that an excluded hidden sector is only made 
more excluded by the addition of heavy degrees of freedom. Note, however, that Fig. [1] also 
shows that in some cases, rather counter-intuitively, an excluded hidden sector may be made 
allowed by the addition of light degrees of freedom. 
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FIG. 1: Bounds from BBN in the (g^ N , g^y) P^ne, where g^ N and g£ e ™ N are the 
hidden degrees of freedom with masses m < T BBN and T BBN < m < T BH , respectively, for 
£rh = Jrjj/Trh = 0.5, 0.7, 0.8, 1.0 (from top to bottom). The regions above the contours 
are excluded. We assume that the observable sector reheats to a temperature above the mass of all 
MSSM particles. The values of (g^^ , <7heavy0 are marked for four example hidden sectors: (A) 
1-generation and (B) 3-generation flavor- free versions of the MSSM with T BBN < mx < ?rh' an< ^ 
(C) 1-generation and (D) 3-generation flavor-free versions of the MSSM with m x < T$ m /2 (see 
Sec. HID. 



CMB temperature anisotropy measurements provide another constraint on the number 
of relativistic species. WMAP 5 year data, combined with distance information from baryon 
acoustic oscillations, supernovae (SNIa) and Hubble constant measurements, imply N e g = 
4.4±1.5 (68% CL) In contrast to the BBN bound, this is a 68% CL bound; at 95% CL, 
the bound is very weak. This bound is therefore numerically much weaker than the BBN 



bound. In the future, however, CMB data from, for example, the Planck satellite [40| will 



achieve a much higher sensitivity in constraining the amount of radiation energy 
441 ]. At the same time, of course, the CMB bound constrains the energy density at much 
later times with redshift z ~ 1000 — 3000 and observable sector temperatures T ~ 1 eV. It 
is thus interesting to explore the CMB constraint when the hidden sector includes particles 



with mass between 1 MeV and 1 eV. (See also Ref. [45J for a combined constraint.) 
The current CMB constraint in our framework is 

Mmb) = I ■ 2 • (iVeff - 3.046) < 1.30 (68% CL) , (14) 

where T^/T 7 = (4/11) 1 / 3 is the neutrino to photon temperature ratio after the era of e ± 
annihilation. The fact that the SM neutrinos were not fully decoupled at e ± annihilation 
is taken into account by using the modified number of SM neutrinos of 3.046 instead of 
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FIG. 2: As in Fig.[TJ but for bounds from the CMB in the (gftgf B , 5^ B ) plane, where gfrgf 5 and 
(/heavy 3 are the hidden effective degrees of freedom with masses m < T^ MB and T^ MB < m < T BH , 
respectively. The values of (g^S^, <?heav^ B ) are given for four example hidden sectors: (C) 1- 
generation and (D') 3-generation flavor-free versions of the MSSM with Tq MB < nix < ^rh> ano ^ 
(C") 1-generation and (D") 3-generation flavor-free versions of the MSSM with mx < Tq MB /2 (see 
Sec. mj). 



3 46|, |47J . Following the BBN analysis, this CMB bound implies 

3 

/ ;, ( \ 3 ,, , ., 
flight 



h CMB < 



4.14 



Heavy 



RH 



/ /iCMB\t frCMB 
flight ) 



(15) 



where g^™ B and gfc^ B are denned as in Eq. ([IT]), but with T^ BN replaced by T£ MB . 

We show current CMB bounds in the (g^™ 3 , S'heavy 6 ) plane in Fig. [2, along with four 
example models. Figure [TJs Model C, with all degrees of freedom light compared to Tg BN , 
is now differentiated into many models, depending on how many of the degrees of freedom 
are light compared to Tq MB . These define a line in Fig. [21 with endpoints given by Models 
C and C". We find that if these models are allowed by the BBN constraints, they are also 
allowed by the current CMB constraint. Similarly, Fig. [TJs Model D is differentiated into 
Models D' and D" in Fig. [2j but again we find that in such models, if the BBN constraint is 
satisfied, the current CMB constraint is also. 

In these models, then, the current CMB constraint is never as stringent as the BBN 
constraint. This may change with the strengthening of the CMB constraint in the future. In 
addition, it is interesting that in the presence of hidden sectors, there may be many particles 
with mass between T BBN and T^ MB , and BBN and the CMB yield independent information 
about the hidden sector mass spectrum. This contrasts with the case of the SM, where CMB 
constraints may be viewed as consistency checks on BBN results. 
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III. A CONCRETE MODEL 



In Sec. [111 we found that hidden sectors may contain many degrees of freedom without 
violating cosmological bounds. With this as motivation, in the following section, we will 
examine freezeout in hidden sectors whose temperatures may differ from the observable 
sector's. To prepare the way, in this section we first specify a model and mass spectrum for 
the hidden sector. 



The WIMPless models of Ref. [30| are supersymmetric, with supersymmetry breaking 
mediated by gauge interactions. A worthwhile endeavor would be to build a complete 
GMSB WIMPless model. In this study, however, we instead work with a "GMSB-inspired" 
model. We do this for two reasons. First, our primary aim here is to clarify cosmological 
issues. Building a complete GMSB model with a WIMPless hidden sector, complete with 
\i term solution and other particle physics details, would take us too far afield, especially 
since, as we will see, the relic density is independent of many of these details. Second, 
as noted in Sec. HJ the essential virtue of WIMPless models, that they have the naturally 
correct thermal relic density, relies solely on the fact that masses in the hidden sector are 
proportional to gauge couplings squared. This is true for GMSB, of course, but is also 
valid in many other contexts, such as anomaly-mediated supersymmetry breaking models, 1 
lessening the motivation for detailed consideration of any one GMSB model. 

The model we consider has an MSSM-like hidden sector, with SU(3) x SU(2) x U(l) gauge 
interactions. In contrast to the MSSM, however, we assume that all Yukawa couplings are 
(9(1). We refer to this as a flavor-free version of the MSSM. For simplicity, we will focus on 
a 1-generation flavor-free version of the MSSM, but, as all Yukawa couplings are large, we 
will use 3rd generation nomenclature for this one generation of hidden matter fermions and 
sfermions. 

All hidden superpartners are assumed to have mass ~ mx, the superpartner mass scale. 
In addition, since nix also sets the scale for the Higgs boson parameters, the hidden Higgs 
bosons h h , H h , A h , and H h± also have this mass, as do the hidden weak gauge bosons W h 
and Z h . Given 0(1) Yukawa couplings, the hidden fermions r h , b h , and t h are therefore 
also at the superpartner mass scale. The light, effectively massless particles are the hidden 
photon / y h , gluon g h and neutrino v h . As is typical of GMSB, the gravitino is also light, 
and its relic density from late decays to the gravitino is negligible. We further assume low 
or moderate reheating temperatures so that the abundance of gravitinos produced during 
reheating is also insignificant. 

WIMPless dark matter requires that one particle at the mx mass scale be stable. We 
assume that the lightest mx scale particle is the hidden right-handed stau X = f^. This 
is the lightest superpartner in much of GMSB parameter space, and it is also reasonable to 
assume that it is lighter than the r h , given 0(1) Yukawa couplings, just as the stop may be 
lighter than the top in the MSSM. The stau's stability is then guaranteed by hidden sector 
electric charge conservation. One might worry about the viability of a charged dark matter 
candidate. Efficient Compton scattering of hidden staus on hidden photons Tr ± 7 /i — » t^ ± 7 /i 
would lead to dissipative energy loss, making the dark matter particles behave more like 
baryons. However, this scattering cross section is cr ~ g hA /m x , which, given Eq. ([2]), is as 
weak as the SM weak interactions. Hidden sector charged dark matter is therefore consistent 
with the known properties of dark matter from structure formation. 



We thank Raman Sundrum for making this point. 
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Given the few light states, there are only a few dark matter annihilation channels: 

f£ + f£ _ -> z/V, 7 V\ 7^ • (16) 

The neutrino channel is mediated by an s-channel Z h . The gauge boson final states are 
produced through processes with fjj in t- and w-channels and 4-point ff77 and ffjZ ver- 
tices. The relic density is therefore dependent on only a few of the many supersymmetric 
parameters, which we will now discuss. Note that we assume no left-right stau mixing. 

All three gauge couplings enter the analysis. The SU(2) and U(l) gauge couplings enter 
the annihilation cross sections. The SU(3) gauge coupling must also be specified, because 
it determines the number of relativistic degrees of freedom at BBN. For simplicity, for the 
hidden gauge couplings, we adopt assumptions motivated by unification. Gauge coupling 



unification implies that the hidden weak mixing angle is tan 6^ = y 3/5 at the GUT scale. 
In principle, this is modified by renormalization group (RG) evolution. The RG equations 
for the hidden sector MSSM gauge couplings are 

^ = J—b h q h3 (17) 

d\n(Q/M) 16vr 2 a Ua ' y } 

where b h a = (13/5,-3,-7) for the 1-generation case, and b h a = (33/5,1,-3) for the 3- 
generation case. As we will see in Sec. IIV} however, WIMPless models are valid in the 
ranges keV ^ mx < TeV and 1CT 5 ^5 gx{mx) < 1, depending on the reheating temperature 
ratio £rh- F° r most of this range, then, gx( m x) ~ 0.1, and the fractional change in the 
gauge couplings from RG evolution is \Ag%/g%\ « \b^\(g^ 2 /16tc 2 ) \n(M mcss /mx) ~ 0.01; that 



is, the gauge couplings evolve very little. Given this, we will assume tan 6^ = y3/5 at m X - 
For mx ~ 100 GeV, this is not the prediction of unification, and our assumption therefore 
implicitly assumes some more complicated structure. 

The unification assumption also typically implies that the hidden confinement scale 



A^ CD ~ mx e b l4 2 (™x) (1 8 ) 

is far below Tg BN . For example, for the 1-generation case, taking mx ~ 100 GeV and 
g^imx) = 0.5, we find Aq CD ~ 10" 9 eV. We therefore assume Aq CD < Tg BN , and so the 
hidden gluons contribute to g^ at BBN. Again, this does not follow from unification for the 
largest m x . 

We will identify gx with the hidden SU(2) gauge coupling g h . As for the hidden MSSM 
mass spectrum, we assume that all massive hidden particles have masses in the range mx < 
rrii < 2mx, and that all relativistic degrees of freedom are always at the same temperature. 
The effective relativistic degrees of freedom in the hidden sector are then 

h h _ h h _ f 116.25 (228.75) , T h > 2m x nq , 
9AJ- )-9*s{-L )-| 19.75(23.25), T h < m x , {J) 

for the 1-generation (3-generation) case, and we linearly interpolate between these two values 
for T h between mx and 2m x - The exact spectrum in this range is largely irrelevant for the 
relic abundance calculation, since we consider only dark matter that is non-relativistic at 
freezeout. Given these assumptions, the thermal relic density is completely determined by 
only four parameters: 

rph 

m x , gx, m z , £rh = 7p~~ ■ (20) 
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T [GeV] 

FIG. 3: Evolution of £ = T h /T, the ratio of hidden to observable temperatures, for mx = 1 GeV, 
1 MeV, 1 keV and 1 eV, assuming Trh = 50 TeV, £rh = 0.8, and a hidden sector that is the 
1-generation flavor-free version of the MSSM. For other £rh> these curves are simply re-scaled by 
£rh/0.8. All supersymmetric particles in the observable sector are assumed to have mass around 
1 TeV. 



In Fig. [3] we plot the relative evolution of the observable and hidden photon temperatures 
for the cases mx = 1 GeV, 1 MeV, 1 keV and 1 eV. To make this plot, we must determine 
T h as a function of T for a given mx- To get the hidden sector temperature T h at a given 
observable sector temperature T, the value of g*s(T h ) is required, which is in fact unattain- 
able unless one knows T h . Our approach is the following: we first solve for g^ s (T h ) T h3 from 
Eq. (|HJ) with T B bn and Tg BN replaced by T and T h , respectively. We then determine T h 
by dividing this quantity by g*s(T h ) as determined by Eq. (TH?]) in the three trial intervals 
T h < mx, mx < T h < 2mx, and 2mx < T h , and choose the value of T h that yields a 
self-consistent solution. 

The general feature of Fig. is that as T decreases from T RH , £ = T h /T drops when 
observable degrees of freedom become non-relativistic and rises when hidden degrees of 
freedom become non-relativistic. Rapid drops occur at T ~ TeV when the observable 

superpartners become non-relativistic and £ oc ^g*s(T), and at T ~ 250 MeV, where the 
QCD transition occurs. The rise occurs when the hidden mx mass particles become non- 
relativistic at T ~ 4mjf — 10mx- At T ~ mx — 5mx, where the exact value depends on 
mx, the hidden sector phase transition concludes, and £ begins to track vV g^siT) again. 



IV. THERMAL RELIC DENSITIES 

We are now ready to determine thermal relic densities in the case of observable and hidden 
sectors with unequal temperatures. For the hidden sector, we assume the 1-generation flavor- 
free version of the MSSM described in Sec. IIHI The temperature of the hidden sector thermal 
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bath at any instant can be read off from Fig. [31 

The number density of dark matter particles n is determined by the competition be- 
tween the expansion rate of the universe and the dark matter annihilation rate through the 
Boltzmann equation 



(IT) 

— = -3Hn-(a A v)(T h ) 
at 



n 



<iT h ) 



(21) 



The Hubble parameter depends on the effective relativistic degrees of freedom in both the 
observable and hidden sectors and is 



H(T) 



4tt 3 G 



N tot 



45 



gT\T)r 



(22) 



where 



(23) 



In contrast to H, (<Jav), the thermally-averaged product of the total annihilation cross 
section and the M0ller velocity, and n eq , the equilibrium number density, are determined by 
the hidden temperature alone. 

To solve the Boltzmann equation, it is standard to change to dimensionless variables 
through 

(24) 



t 



x 



T 



n 



Y 



n 
s 



These new variables are natural from the observable sector viewpoint, as they use the ob- 
servable sector's temperature T as the "clock," and the observable sector's entropy density 
s = %- g*s T 3 as the fiducial quantity characterizing the universe's expansion. Note, how- 
ever, that in the present context there is also another useful change of variables, 



t 



h _ m x 

X - J, h 



n 



Y h = - 
s h 



(25) 



that is far more natural from the hidden sector point of view. We will return to this below. 
In terms of the observable sector variables of Eq. ( 1241) , the Boltzmann equation becomes 



dY -A 

dx 



Y 2 



with the expansion and annihilation effects encapsulated in the quantity 

9*s(T) ( 1 T dg* s (T)\ m x 



A 



7T 



45 G 



1 



3 0* S (T) dT 



{a A v)(T h 



(26) 



(27) 



We will solve this equation in the non- relativistic regime, where x h > 3. The equilibrium 
value of the number of dark matter particles in a comoving volume is 



Y eq (x h ) = 



n eq {T h 



1 9 



s(T) " s(T) 2tt 2 J mx e E / Th - 1 



E 2 — m\ 



EdE 



45x 



gK 2 {x h )e 



4tt 4 ^ s (T) 



(28) 
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where g denotes the number of internal degrees of freedom of the dark matter particle, and 
the last approximation of Eq. (1281) is valid for x h ^ 3. Following Ref. [50], the thermally- 
averaged cross section times velocity is determined using 



8m x T h K$(^) 



jT dsa^s^is-Am^K^] . (29) 



In Eqs. (1281) and (129|) . Ki are modified Bessel functions of order i. Equation (129!) was 
derived for particles in thermal equilibrium obeying Maxwell-Boltzmann statistics, but it is 
applicable to other statistics, provided T h < m x /3. 

We determine the fjj total annihilation cross section cta(s) with the help of the CalcHEP 
package (H, 49[. Expanding around e = (s — 4m x )/4m x (see Ref. 50]), we find 



(a A v)(T h 



rn 



x 



a + ai 



( 



m x 



V T h 



+ a 2 



\ T h) 



+ ... 



(30) 



where 



a i 




3 
2 



771 



4m 2 



tan 2 e% 



1 



-— — tan 2 6^ + 

6vr 3vr w 12vr 



sin 4 ^ 



cos 4 ^ 



-4 + 



z h 



+ 



(31) 



sin 4 ^ , 



(32) 



and r Z h is the Z h boson decay width. In keeping with our non-relativistic approximation, 
we neglect the a 2 and higher order terms. For m Z h w 2mx, t l r + t l r~ v h v h is greatly 



enhanced by the s-channel resonance. Off this resonance, however, the neutrino channel 
is typically sub-dominant, because it is P-wave suppressed, as can be seen above. In our 
numerical study, we fix m Z h = 1.5 mj. The Z h boson remains in equilibrium with the 
hidden thermal bath through the decay and inverse decay processes Z h <-> v h v h . 

We now present several results for Y(x). Because of the "stiffness" of the Boltzmann 
equation of Eq. ([26]) . we adopt the implicit trapezoidal method with adaptive step size, 
developed by the authors of the DarkSUSY package [5l[, to solve for Y(x) numerically. We 
first examine the £rh dependence of the thermal relic density. In Fig. HI we plot Y(x) for 
various values of £rh, m x = 1 MeV, and g x = 5.57 • 10~ 4 . The latter two parameters 
are naturally accommodated in the WIMPless GMSB framework and yield Q x h 2 in the 
required range, as we will see below. As £rh decreases, we see two effects. First, for a given 
x, the hidden sector is colder, and so Y eq (x) is smaller. Second, since freezeout occurs when 
H ~ T = n(aAv), when £ RH decreases and n drops, freezeout occurs earlier. As can be 
seen from Fig. HI these effects act against each other, and the thermal relic density is fairly 
insensitive to £rh- 

These results may be understood more precisely using standard analytic approximation 
formulas for thermal relic freezeout. The equilibrium number density is 



Y 



eq 



O.U5{g/g* s )x 3/2 e /2 e 



(33) 
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(m x =1 MeV, g x =5.57e-4) 



| RH =0.8 
| RH =0.5 
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x = m x / T 
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FIG. 4: Y(x), the fjj number density per comoving volume, as a function of x = mx/T, for 
mx = 1 MeV, gx = 5.57 • 10~ 4 , and various ratios of hidden sector to visible sector reheating 
temperature £rh = 0.8, 0.5 and 0.3 (from top to bottom). 



where £ = Th/T, as defined above. Following the procedure given in Ref. and assuming 
5*-wave annihilation, an approximate expression for x at freezeout is 



Xf ~ £ In 



0.038 M P1 m x a (g/Jg^) £ 3/2 6(6 + 2) 



0.038 M P1 m x a {g/Jg«*) £ 3/2 6(6 + 2 



(34) 



where o"o = aog x /m x and the parameter 6 is tuned to make these analytical results fit the 
numerical results. 2 

For iS-wave annihilation, the final relic abundance is approximately 



3.79 a;/ 



9*s/^fg¥ t ) M PI m x a 



(35) 



which depends linearly on £ through Xf. This linear dependence is verified by the numerical 
results presented in Fig.|H The relative values of Y obtained from Eq. (1531) are 3.6 : 1.9 : 1 for 
£rh = 0.8, 0.5, 0.3, respectively, which almost exactly match the results from our numerical 
computation, which yields 3.5 : 1.9 : 1. 

The dependence on £rh may also be understood from the hidden sector viewpoint. In 
Fig. [5l we consider the same model parameters as plotted in Fig. HI but we now present 
Y h (x h ). From a hidden observer's point of view, the dependence on £rh arises solely because 
£rh determines how hot the observable sector is, and this impacts the expansion rate H. For 



2 For £rh = 0.3 (0.8), we find that choosing S = 0.2 (0.5) yields agreement typically better than 3%. For 
larger mx or higher £rh, the final relic density is less sensitive to S. 
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FIG. 5: As in Fig. 01 but now plotted in terms of the hidden sector parameters Y h (x h ). 



a fixed x , lower £rh implies higher T, faster expansion, and earlier freezeout. This is seen 
in Fig. [5j At first sight, the result that lower £rh implies larger Y h may seem to contradict 
the results obtained earlier, where lower ^rh implies lower relic density. Note, however, that 
to convert Y h to a relic density, one must multiply Y h by the current s h , which, unlike s, 
is not independent of £rh- For lower £ R h, the current hidden sector temperature is lower, 
and the current s h is also lower. This effect makes the relic density smaller for lower £rh, in 
accord with our previous results derived from the observable sector viewpoint. 

We now examine the dependence of the thermal freezeout process on mx and gx- In 
Fig. El we plot Y(x) for fixed £rh, but for a wide range of (mx,gx) with mx oc g x , the 
scaling relation of WIMPless scenarios. For low mx, freezeout occurs earlier, as expected 
given Eq. (CM)) . In our analysis, we have assumed that freezeout occurs in the non-relativistic 
regime, where x h > 3. Converting x h to x using Fig. [21 we see from Fig. Elthat our analysis 
is self-consistent even for mx as low as 10 keV. 

In Fig. [7] we present contours of constant Qxh 2 in the (mx,gx) plane. These contours 
result from our numerical analysis; although we have confined our previous plots to x ^ 100, 
the results of Fig. [7J use the relic density as determined at x ~ 1000, where it has truly 
stabilized. 

Figure [7J summarizes all of our work so far, and there are several noteworthy features. 
First, we see that these contours essentially follow the scaling relation mx oc g\. If, for 
example, mx is lowered by a factor of 4 and gx is lowered by a factor of 2, the weakened 
interaction implies earlier freezeout and larger Y Q , but this compensates the lower mx to keep 
the relic density Qxh 2 oc mxY approximately invariant. Alternatively, using the analytic 
results presented above, this scaling relation keeps a a invariant, and, given Eq. (135]) . this 
also keeps Qxh 2 invariant. Of course, this is as expected from the general arguments given 
in Sec. HI but Fig. [7J shows that these arguments, based essentially on dimensional analysis, 
are in fact supported by a detailed calculation in the context of a concrete model with dark 
matter masses varying greatly and far beyond the conventional WIMP range. 
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m x =1 GeV, g x =0.0176 
m x =1 MeV, g x =5.57e-4 
m x =10 keV,g x =5.57e-5 




(^RH=0-5) 



>- 



10 



J L 



1 



10 



100 



x=m x / T 



FIG. 6: Y(x), the fjj number density per comoving volume, as a function of x = mx/T, for fixed 
£ R H = 0.5 and (m x ,gx) = (1 GeV, 1.76 • 1(T 2 ), (1 MeV, 5.57 • 1(T 4 ), and (10 keV, 5.57 • 10~ 5 ). All 
cases are chosen to have the same ratio rax I g\ ■ 

Second, for comparison, Fig. [7] also includes two dashed contours to indicate the corre- 
sponding observable sector weak mass scale m wea k- The definition of m wea k depends on how 
one chooses to characterize the natural weak scale in supersymmetric theories (the \x param- 
eter, the Bino mass, etc.), and these depend on several parameters, such as the number of 
messengers and the messenger scale. We simply choose m weak = (m x /g x )g' 2 , where g' is the 
U(l) gauge coupling of the MSSM, and plot contours for m weak = 100 GeV and 1 TeV. The 
fact that there exist values in the (mx,gx) plane that give the correct flxh 2 is, of course, 
not surprising. However, the fact that the two solid curves lie between the two dashed curves 
means that the preferred values of (m x ,9x) are those that correspond to models in which 
the observable sector's superpartners are at the weak scale. 

We also see from Fig. [7] that, for a fixed m x , lower £ R h requires lower g x to obtain the 
correct relic density. This is because, at a given observable temperature T, colder hidden 
sectors imply lower equilibrium number densities. To keep the relic density constant, gx 
must be lower so that freezeout happens earlier. This is consistent with the results shown 
in Fig. |H 

Finally, note that the £rh = 0.3 contour ends at mx — 45 keV. This is because the 
criterion that dark matter freezeout while non-relativistic, that is, x j ^ 3, is violated below 
this mass. In Fig. [BJ we plot, as a function of £rh, the lower mass limit resulting from 
requiring that freezeout occurs with x h j ^ 3 and yields Qh 2 = 0.11. We can see that the 
lower mass limit goes up with smaller £rh- This is because for colder hidden sectors, freezeout 
occurs earlier, and to satisfy x*j > 3, we need larger m X - This exercise establishes that the 
WIMPless dark matter framework may be valid at least down to dark matter masses of 
mx ~ keV. Below this mass, freezeout is not non-relativistic, and our analysis breaks down. 
It would be interesting to perform a more precise analysis without assuming non-relativistic 
freezeout, but this is beyond the scope of our study. 
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FIG. 7: Contours of £lxh 2 = 0.11 in the (mx,gx) plane for £rh = 0.8 (upper solid) and 0.3 (lower 
solid), where the hidden sector is a 1-generation flavor- free version of the MSSM, corresponding to 
Model A (for mx > 1 MeV) or Model C (for mx < 1 MeV). The minimum mx consistent with 
the requirement of non-relativistic freezeout depends strongly on £rh (see text). Also plotted are 
lines of m wea k = i m x / gx)g' 2 = 100 GeV (upper dashed) and 1 TeV (lower dashed). 

V. IMPLICATIONS OF CONNECTORS 

Throughout this work, we have ignored the possible existence of connector particles Y 
with both SM and hidden gauge charges. If these are added, do they significantly affect the 
conclusions derived above? 

There are two significant roles that connector particles can have. First, they likely estab- 
lish thermal connectivity between the observable and hidden sectors when the temperatures 
are above the connector mass scale my. At temperatures below my, however, the hidden- 
SM interactions become weak, and the observable and hidden sectors decouple. As shown 
in Figs. [T] and [2J the BBN and CMB constraints actually allow, for example, a 1-generation 
flavor-free version of the MSSM with mx ^bbn; even when £rh = 1- In this case, then, 
there is no problem if the two sectors are in thermal contact right after reheating. For other 
hidden sectors, it may be possible that the existence of connectors enforces £rh = 1> or & lter- 
natively, the assumption of £rh 7^ 1 implies an upper bound on the reheating temperatures. 
One must check on a case by case basis whether connector particles upset the premises of 
our calculations. 

A second affect is that connector particles open new avenues for dark matter annihilation. 
As an example, we assume that fermionic connector particles Y couple through Yukawa 
couplings XXYf, where / is a SM fermion. We expect m Y ~ max(mj,m weak ) [30(], and also 
require m Y > m x to prevent the decay X — > Yf. 

Such connector interactions induce annihilations XX — > // through processes with Y 
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FIG. 8: Lower bound on the dark matter particle mass mj for different £rh in the WIMP less 
scenario. As in Fig. [7J the hidden sector is a 1-generation flavor-free version of the MSSM, corre- 
sponding to Model A (for m x > 1 MeV) or Model C (for m x < 1 MeV). 



particles in the t-channel. There are two cases. If m x ^ ^weak, m Y ~ ^weak, and 

. , \ A m\ , A 4 A 4 
m Y m Y m 2 weak 

On the other hand, if m x ^ ^ W eak, m Y ~ m x, and 

A 4 A 4 

(<Tv)f f r < — — . (37) 

m x m4 eak 

In either case, the annihilation cross section to SM particles is small relative to the completely 
hidden annihilation cross sections ^4eak/ m weak> provided A ^ g WC ak — 0.65. This is a rather 
weak criterion, and so typically, annihilation to SM particles is sub-dominant, and the 
analysis given above holds even in the presence of connector particles. 



VI. CONCLUSIONS 

In the paper, we have studied the possibility that dark matter is hidden, that is, has 
no SM gauge interactions. Hidden sectors appear in many frameworks for physics beyond 
the SM, and hidden dark matter is perfectly viable, given all observations to date. We are 
particularly motivated by the recent proposal of WIMPless dark matter. As with WIMPs, 
WIMPless dark matter has the key virtue that its thermal relic density is in the right range 
to be cold dark matter. In contrast to WIMPs, however, WIMPless dark matter may have 
masses and couplings that differ drastically from WIMPs. 

Although hidden sectors interact very weakly with the SM, hidden particles contribute 
to the energy density of the Universe. This impacts the expansion rate of the Universe, and 
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so bounds from BBN and the CMB constrain hidden sectors. As is well known, current 
constraints completely exclude a hidden sector that is an exact copy of the MSSM. We find, 
however, that these constraints are rather brittle and easily avoided if the hidden sector is 
just slightly colder than the observable sector or if the mass scale of hidden superpartners 
differs from the MSSM. Several viable examples are discussed in Sec. [TT1 and Figs. [T] and [2] 
give model-independent bounds that can be used to determine the viability of other models. 

In Sec. 1111} we then defined a concrete model, a 1-generation flavor-free version of the 
MSSM. This model has features generic to GMSB models and provides a WIMPless dark 
matter candidate, the hidden stau fjj. The fjj relic density depends sensitively on only a 
small number of parameters. This relic density was examined in Sec. II VI by solving the Boltz- 
mann equation numerically, and understanding these results through well-known analytic 
approximations, generalized to hidden sectors with different temperatures. These results 
confirm that in WIMPless models with mj oc g x , the relic density is naturally in the right 
range for dark matter masses keV < mx ^ TeV, greatly extending the conventional WIMP 
mass range. WIMPless dark matter therefore provides a class of dark matter candidates 
that shares the key relic density virtue of WIMPs. This generalization may have interesting 
applications, given the diverse and interesting phenomenology and observational anomalies 
already present at the keV, MeV, GeV, and TeV mass scales. 

We have determined the allowed range of mx by requiring that this dark matter freezes 
out with the right relic density while non-relativistic. Of course, this range may be further 
constrained by other considerations. In particular, the lightest WIMPless candidates, with 
mx ~ keV, may have interesting implications for structure formation. After freeze out, 
they can still couple to the hidden sector thermal bath through elastic scattering processes 
T~Ftf h ~ * Tftf and Tfti> h — > f\v h . After their kinetic decoupling at temperature T^, they 
then stream freely between over- and under-dense regions. Free-streaming of thermal relics 
causes damping in density perturbations (52| at a comoving scale Afs oc m x X ^ T v ^ 2 \ in 



contrast to WIMPs [531. 1541] . this is the dominant effect for MeV or lighter dark matter [55 
If this scale falls in the range (1 — 80) h^ 1 Mpc, the linear matter power spectrum from 
high-resolution Lyman-a forest data can be used to place a lower bound on the mass of the 
dark matter particle. Current lower mass limits on warm dark matter of 0.55 — 0(1) keV 
from Lyman-a data [56], 57, 58, 5^] therefore imply lower mass bounds in the WIMPless 



scenario, which depend on the kinetic decoupling temperature T^d for different £rh- 

In addition, another kind of dark matter lower mass limit follows from the observation 
that the phase-space density of the dark matter particles in galaxy halos cannot exceed the 
maximum value of the phase-space density when the dark matter particles were in kinetic 
equilibrium [60]. With knowledge of the WIMPless dark matter's chemical potential, one 
can then apply these generalized Tremaine-Gunn bounds from Ref. [6lJ to our scenario. 
Detailed considerations of these effects are beyond the scope of this study, but they may 
modify the lower bound on mx or, alternatively, provide interesting astrophysical signals 
for the hidden dark matter scenarios discussed here. 
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